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Abstract 

We reinvestigate the QCD sum rule for the ttNN coupling constant, g, start- 
ing from the vacuum-to-pion matrix element of the correlation function of the 
interpolating fields of two nucleons. We study in detail the physical content of 
the correlation function without referring to the effective theory. We consider 
the invariant correlation functions by splitting the correlation function into 
different Dirac structures. We show that the coefficients of the double-pole 
terms are proportional to g but that the coefficients of the single-pole terms 
are not determined by g. In the chiral limit the single-pole terms as well as 
the continuum terms are ill defined in the dispersion integral. Therefore, the 
use of naive QCD sum rules obtained from the invariant correlation functions 
is not justified. A possible procedure to avoid this difficulty is discussed. 
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The QCD sum rule invented by Shifman, Vainshtein and Zakharov provides us with a 
way to relate the physical quantities of the hadrons to the matrix elements of the quark 
gluon composite operators by means of the operator product expansion (OPE) 0,3- 

There has already been quite a history in the study of the meson-baryon-baryon coupling 
constants within the framework of the QCD sum rule 0-0. Two ways of the formulation 
exist in constructing sum rules: one is to start from the vacuum-to-vacuum matrix element 
of the correlation function of the interpolating fields of two baryons and one meson and 
the other is to start from the vacuum-to-pion matrix element of the correlation function 
of the interpolating fields of two baryons. In the pionnering work by Reinders, Rubinstein 
and Yazaki the former approach was taken but in the following works including the one by 
the same authors the latter approach was employed. This is because the former approach 
has to assume additional extrapolation from the large space-like momentum to the zero 
momentum for the pion compared to the latter approach. Furthermore it was shown in 
Ref . that the three-point function method for the treatment of both the isospin conserving 
and isospin violating pion-nucleon-nucleon couplings is plagued by problems with higher 
resonance contamination. 

Therefore, we only discuss here the formulation with the vacuum-to-pion matrix element 
of the correlation function of interpolating fields of two baryons: 

n(p,Jfe) =-ij d 4 xe ipx (0\T{ri(x)f)(0)]\n(k)), (1) 

where T}{x) is the interpolating field for the nucleon, |7r(fc)) is the pion state with momentum 
k normalized as {ir(k')\7r(k)) = 2k (2ir) 3 5 3 (k' — k) and the isospin is neglected for simplicity. 
The correlation function can be split into four Dirac structures as, 

1% k) = ij^ip^pk) + ll5 ^U 2 (p 2 ,pk) + ll5 tu 3 (p 2 ,pk) - ij 5 m ^ k U 4 (p 2 ,pk), (2) 

mm m z 

where m denotes the nucleon mass, IT (i = 1 ~ 4) are functions of Lorentz invariants, p 2 
and pk, and are called invariant correlation functions. Reinders, Rubinstein and Yazaki 
constructed a sum rule for the pion-nucleon-nucleon (nNN) coupling constant in the chiral 
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limit with k = from the invariant correlation function with the Dirac structure 275, III, 
in the leading order of the OPE. They showed that if the sum rule is divided by the odd- 
dimensional sum rule for the nucleon mass, the result is consistent with the Goldberger- 
Treiman relation with gA = 1, where gA is the axial charge of the nucleon 0. Shiomi 
and Hatsuda improved the sum rule by taking into account higher dimensional terms and 
a s corrections of the OPE and also the continuum contributions [|J. They showed that 
even after these corrections are taken into account, the Goldberger-Treiman relation with 
gA = 1 holds as long as the same continuum threshold is taken in the sum rule for the tiNN 
coupling constant and the odd- dimensional sum rule for the nucleon mass. In Ref. || Birse 
and Krippa pointed out that in the chiral limit the vacuum-to-pion correlation function is 
obtained just by chirally rotating the vacuum-to- vacuum correlation function and therefore 
it is obvious that the ratio of the two sum rules is consistent with the Goldberger-Treiman 
relation, which is just a consequence of the chiral symmetry. Then, they tried to obtain a 
new sum rule for the nNN coupling constant which is not just a consequence of the chiral 
symmetry. They considered the invariant correlation function with the structure 275^ and 
obtained a sum rule by taking the limit, k = 0, after removing 275^. In all of these works 
the sum rule is constructed by making an "ansatz "for the absorptive part of the correlation 
function based on the effective lagrangian with the pseudoscalar coupling scheme. In Ref. || 
Kim, Lee and Oka examined how the choices of the effective lagrangian, i.e. the pseudoscalar 
and pseudovector coupling schemes make differences in the sum rules for all the four Dirac 
structures. They concluded that only the invariant correlation function with the Dirac 
structure ^cr^p^k 1 ' , is independent of the two coupling schemes and obtained a sum rule 
by taking the limit, k = 0, after removing j^a^p^k^ . They also claimed that their sum rule 
suffers from less uncertainties due to QCD parameters. 

Even though we agree with Kim et al. in the point that the sum rule should not depend 
on the choice of the coupling scheme, we find their argument unsatisfactory because they still 
rely on particular effective lagrangians. Actually, we can identify the physical structure of 
the correlation function just by general principles without referring to effective lagrangian |§. 
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The purpose of the present paper is to reinvestigate the QCD sum rule for the nNN 
coupling constant by applying the procedure of Ref. ||. 

We clarify the physical contents of the correlation function defined by Eq. (|1]) without 
referring to the effective theory. The correlation function II (p, k) in Eq. ([I]) has a pole at 
p 2 = m 2 and (p — k) 2 = m 2 where p and p — k, respectively, become on-shell momenta for 
the nucleon. The -kNN coupling constant, g, is defined through the coefficient of the pole 
as 

u(pr) (p 1 - m)n(p, k) (/} - $ - m)u(qs) | p2=m 2 i(p _ fc)2=rn 2 = -igu{pr)lhu{qs) , (3) 

where q = p — k, u(pr) is a Dirac spinor with momentum p, spin r and is normalized as 
u(pr)u{pr) = 2m and the unnormalized nucleon interpolating field, rj, is replaced by the 
normalized nucleon field, ip. 

In addition to the pole singularity, II (p, k) has a branch cut singularity starting from the 
threshold of the pion-nucleon channel, p 2 = (m + m n ) 2 , to infinity where is the pion 
mass. We restrict ourselves to the positive energy region, p > 0, for simplicity. 

In order to classify the singularity of U(p,k), we define the vertex function, T(p,q',k') 
(p = q' + k'), and the pion-nucleon T-matrix, T(q', k', q, k) (q' + k' = q + k), by 

T(p,q',k') = ($- m )[-i J d A xe ipx (0\T[^(x)^(0)}\7r(k')}](i' - m), 

T{q',k',q,k) = {£'-m)[-i J rf 4 xe Vx (7r(A;')|T[^(x)^(0)]|7r(A;))](^ - m). (4) 

The correlation function is related to the vertex function as 

i> + m i> — h + m 
n(p, k) = — ~T(p,p - k, k)- — -. 

Let us regard Tl(p, k) as a function of the center-of-mass energy, i.e. po in the frame p = 0, 
and consider the discontinuity of the correlation function. We adopt the following notation 
for the dispersive (continuous) part and the absorptive (discontinuous) part, respectively: 

ReF(p) = ^[F(p)\ p o =p o +iri + F(p)\ p0=p o_ iv \, 
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lmF(p) = 2-[F(p)\ p0=p o +irj - F(p)\ p0=p o_ 



■ir\\ ■ 



(5) 



The absorptive part of the correlation function, lmll(p, k), can be written as 
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ImU(p, k) = {j> + m) < Im 



-ReT(p,p — k, k) 



(p 2 — m 2 )((p — k) 2 — m?) 

+Re 71 — — -lmT(p,p-k,k)\(^-jk + m). (6) 

(p z — m z )({p — ky — m z ) J 

The first term represents the pole contribution and the second term represents the continuum 
contribution. 

When the cm. energy, p , is above the threshold of the pion-nuclon channel, m + m„ but 
below that of the next channel, only the pion-nucleon channel contributes in the intermediate 
states and the absorptive part of the vertex function, lmT(p,p — k, k), is given by 

7T t r dSV 



Imr(p, p — k,k) 



4p J (2tt)3 



T(p,q',k')(i' + m)T*(q',k',q,k), 



(7) 



where t is defined through po = \J m 2 + t 2 + \Jm 2 + t 2 . q' and k! are the on-shell momenta 
for the nucleon and the pion, respectively. 

From Eq. (0), one sees that in the vicinity of the pion-nucleon threshold T can be ex- 
panded as 



lmT(p,p — k, k) 







(po < m + m n ) 



ReT(p,p — k, k) 



(8) 



Gif + 0(t 3 ) (p >m + m n ) 
G + G x t + 0(r 3 ) (p <m + m 7r ) 
G + O(t 2 ) (p >m + m n ) 

where r = it. In Eq. (§) the behaviour of the dispersive part is determined from that of the 
absorptive part by the Cauchy-Riemann equations. From Eq. (|) one sees that 



_d_ 

dp 

_d_ 

dp 



ReT(p,p — k, k) 
ReT(p,p — k, k) 



p Q =m+m 7T -e 



CO 

finite constant 



(m* + 0) 
(m n = 0) 



(9) 



po=m+m^+e 

which is crucial in the following discussions. 



We discuss the physical content of the invariant correlation function. Similarly to Eq. (0) 
we split the vertex function into four Dirac structures as 

T(p,p - k, k) = i^ 5 Ti(p 2 ,pk) + ij 5 —T 2 (p 2 ,pk) + i^—Y^.pk) - i ^-^K - T±(p 2 ,pk), 

mm m l 

where 1^ ~ T 4 are the invariant vertex functions. 

The invariant correlation functions are related to the invariant vertex functions as 

f (ri - r 2 - 2r 3 + (2 - r 4 ) \ (r, + r 2 - £r 4 ) + ^( P 2 - m 2 )r 4 



ni 



(p 2 — m 2 )((p — k) 2 — m 2 ) (p — k) 2 — m 2 

p2 _ m 2 2m 2 ' 

■j-j- _ r 3 — r 4 ^ r 2 + r 3 — r 4 

2 p 2 — m 2 (p — k) 2 — m 2 ' 

" 2 .^ V.\ ri 1 -n 3i-i 1-n 



n m 2 {-T 1 + T 2 + 2T 3 -{2-^)T,) p 2 + T 3 - |r 4 ±r 

J-J-3 — KTT, 7TT; 7Ts r 



(p 2 — m 2 )((p — k) 2 — m 2 ) (p — k) 2 — m 2 p 2 — m 2 ' 

n ^ ^(r 1 -r 2 -2r 3 + (2-^)r 4 ) | ip 4 | ip 4 _ 

(p 2 — m 2 )((p — k) 2 — m 2 ) (p — k) 2 — m 2 p 2 — m 2 

Noting that the ttNN coupling constant, g, is related to the invariant vertex functions as 

,2 \ 



(171 



. (11) 

p 2 =m 2 , (p— k) 2 =m 2 



we clearly see that the coefficients of the double poles at p 2 = m 2 and (p — k) 2 = m 2 
are proportional to g, i.e. —^f-g, 0, m 2 g and —m 2 g for the structures 275, 275^, ijsjb 
and j5<J fJjU p tJ, k 1 ' ', respectively, but that the coefficients of the single poles at p 2 = m 2 or 
(p — k) 2 = m 2 are not determined by g. 

In the chiral limit with k — the invariant correlation functions can be regarded as 
functions of p 2 and the discontinuities of the invariant correlation functions across the real 
p 2 axis are given by 

p 

Imrix = 7rS(p 2 — m 2 )Ti{m 2 ) -\mFAp 2 ), 

p l — m l 
p 

Imn 2 = nS(p 2 - m 2 )r 2 (m 2 ) -lmT 2 (p 2 ) 

p 1 — m l 

Pf 

Imll 3 = nS'ip 2 — m 2 )m 2 go — ix5(p 2 — m 2 )m 2 ReT' A (m 2 ) + — — —m 2 lmT a(p 2 ) , 

\jr — m 1 ) 1 



ImlLi = —ixd'ip 2 — m 2 )m 2 go + irS(p 2 — rn 2 )m 2 ReT' B (m 2 ) — — —m 2 lmT b(p 2 ) , 

\yr — m z y 

where 

t a = -r\ + ^r 2 + (1 + r 3 - 2^r 4 , 
v B = -r 1 + r 2 + 2r 3 - f 1 + ) r 4 , 

and go denotes the nNN coupling constant in the chiral limit. For the structure ^75 the 
double pole vanishes, while for the structures ij^fc and ^a^p^W the double poles survive but 
275^ and r y5<r^ l ,p IJ 'k u themselves vanish. The coefficients of the single poles for the structures 
275^ and ^^CT^vP^k 11 involve the derivatives of the invariant vertex functions with respect to 
p 2 and are therefore ill defined due to Eq. ([5]). The third term of Imn 3 as well as Imn 4 
behaves as 9(p 2 — m 2 ){p 2 — m 2 )^ 1 since the absorptive part of the vertex function behaves 
as p 2 — m 2 due to Eq. (||), whose integral is also ill defined. We should note that the total 
dispersion integral is well defined, but when it is split into two parts as in Eq. ([|) each part 
is ill defined. 

If one employs the effective lagrangian and calculates the vertex function at the tree 
level, one obtains T 1 = —g and T 2 = T 3 = T 4 = for the pseudoscalar coupling scheme, 
while r 3 = I and I\ = T 2 = T 4 = for the pseudovector coupling scheme. If one substitutes 
these invariant vertex functions in Eq. (0), one finds that two coupling schemes give the 
same result only for ImII 4 . This observation made Kim et al. in Ref. claim that only the 
structure ^^cr^ u p^k v is independent of the coupling schemes. (They also claimed that the 
structure has less uncertainty in the OPE.) 

The tree level calculation is meaningful for the double-pole term because the effect of 
the loops can be absorbed by the renormalization of the coupling constant. But the tree- 
level calculation is not meaningful for the single-pole and continuum terms because the 
effect of the loops drastically change the structure. In fact, at the one-loop level the second 
and third terms of Imll 3 and ImII 4 are generated which behave as expected by the general 
consideration above. 

8 



Another point worth while mentioning is the contribution of the excited baryon which 
couples to the nucleon-pion channel. It is certainly true that the effective lagrangian which 
describes such a coupling generates a finite single-pole term. But it is just a part of the 
contribution and it alone does not provide us with the structure expected from the general 
principle. Thus, just adding the tree level contribution of the effective lagrangian for the 
excited baryon is insufficient. 

Therefore, in our opinion a reasonable procedure to make an ansatz for the absorptive 
part of the invariant correlation function to be used in the sum rule would be first to 
transform it, for instance by multiplying by p 2 — m 2 , to a form which does not include 
ill-defined terms and then to parameterize it in accordance with the structure expected by 
the general principle, not by the tree level of the effective lagrangian. 

Another way to construct a well-defined sum rule, a projrcted correlation function ap- 
proach, is proposed in a separate paper [[|. 

In summary, we have reinvestigated the QCD sum rule for the ttNN coupling constant, g, 
starting from the vacuum-to-pion matrix element of the correlation function of the interpo- 
lating fields of two baryons. We have studied in detail the physical content of the correlation 
function without referring to the effective lagrangian. We have considered the invariant cor- 
relation functions by splitting the correlation function into different Dirac structures. We 
have shown that the coefficients of the double poles are proportional to g but that the co- 
efficients of the single poles are not determined by g. In the chiral limit the double-pole 
terms survive only for the Dirac structures ij^jb and 75CT ^p^ k u . For these structures the 
single-pole terms as well as the continuum terms turn out to be ill defined in the dispersion 
integral. Therefore, the use of naive QCD sum rules for these structures in the chiral limit 
is not justified. A possible procedure to avoid this difficulty has been discussed. 
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